Cosmological long wavelength perturbations 



W. Unruh 

Program in Cosmology and Gravity of CIAR 
Dept Physics and Astronomy 

University of B. C. 
Vancouver, Canada V6T 1Z1 



00 

■ Abstract 

This paper presents an exact solution to the long wavelength perturbations for 
Q . the scalar modes and for a scalar field theory with arbitrary potential. Locally 

these modes are coordinate transformations of the homogeneous background 
solutions (although non- locally they are not). These solutions are then used 
to discuss a couple of recent papers in which such perturbations play a role. 

■ Abramo, Brandenberger, and Mukhanov have recently argued that long wave- 
length perturbations have the effect of driving the cosmological constant to 
zero if the higher order perturbation equation are examined. I argue that 
this effect is invisible to any local observer, and thus does not constitute a 
relaxation of the cosmological constant in the normal sense of the term. 

Grishchuk has argued that the standard lore on the strength of the per- 
O { turbations at the end of inflation is wrong. I discuss the disagreement in light 

of the exact long wavelength solutions, and emphasize the importance of the 
initial conditions in resolving the disagreement. 
In the inflationary models for the growth of the universe, the physics of long wavelength 
perturbations (wavelengths which are much longer the the Hubble radius) play an important 
role. It is the behaviour of the perturbations during this time period which determine the 
effect of the perturbations on the present structure of the universe. Recently, a couple of 
papers have discussed these long wavelength perturbations in different contexts. Abramo, 
Brandenberger, and Mukhanov |IJ have argued that higher order corrections to Einstein's 
equations for long wavelength perturbations have the effect of creating a negative cosmo- 
logical constant which reduces the effective cosmological constant for the space-time. They 
also speculate that this could be a mechanism to drive the actual cosmological constant to 
near zero (or in particular to near the current mass density of the universe). On the other 
hand, Grishchuk [||,|| has claimed that the standard calculations for the size of the density 
perturbations at the end of inflation are wrong. Various authors ||[7| have claimed that his 
treatment of the behaviour of the long wavelength perturbations is wrong. 

In this paper, I derive an exact expression for these long wavelength perturbations in 
a system in which the dominant gravitating matter is a scalar field \I/ with an arbitrary 
potential V(^). I note in passing that these solutions therefor apply without much change 
to the evolution with an arbitrary perfect isotropic fluid, under the usual identification that 

p+p = \<r*,^ v (i) 
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I. GAUGE INVARIANCE VS GAUGE FIXING 

Let me first discuss the ABM paper. I was confused by their result that the long wave- 
lengths could renormalise the cosmological constant since the effect is caused by the averaged 
energy momentum tensor of modes whose wavelength is much larger than the Hubble's ra- 
dius at any time of interest. "How could such long wavelength modes affect the cosmological 
constant as seen by an observer who can measure the cosmological constant only over a re- 
gion smaller than his own Hubble radius?" I would have expected modes, whose wavelength 
is much much larger than the Hubble radius, to look like a homogeneous universe to a local 
observer. 

In addition, as they point out, the effective stress energy tensor of lower order pertur- 
bations can be coordinate dependent ( called gauge dependent if one is looking at small 
coordinate transformations). Their approach to this problem was to write the equations in 
terms of "gauge independent" variables. However, as I will argue in the following, there 
is no difference between such a "gauge invariant" approach, and an approach which fixes 
the gauge in some way. Any variable in a gauge fixed formalism is completely equivalent 
to some gauge independent variable, and vice- versa. In particular, this means that, even 
though the effective stress energy tensor is not independent on gauge transformations in 
the gauge invariant approach, it is dependent on which set of gauge invariant variables one 
chooses. In the gauge fixing approach, the same problem arises in that the effective stress 
energy tensor is dependent on which gauge fixing choice one makes. 

Let me review their approach. Instead of writing Einstein's equations in terms of the free 
metric and matter perturbations, Sg^ u , 5^/, they choose a set of vector fields X^ created 
out of the metric and field variables, (i.e., X M = X^^Sgx^S^) ), and defined new gauge 
invariant perturbations 

= e Cx (g(o)^ v + Sg^ u ) - #(0)^ (4) 

where Cx is the Lee derivative with respect to X. A similar definition applies to Sxjj. X is 
chosen so that the barred quantities are invariant under a (first order) gauge transformation. 
Given the free first order metric perturbations written as 

^ \Si + B,i <pg(o)ij + Q\ij + Ri\j + hij 

where 0, B, andQ represent the scalar perturbations, the transverse vector fields,^, Ri rep- 
resent the vector perturbations, and the transverse-traceless tensor represents the gravi- 
tational wave perturbations, g^)^, are assumed to be a background flat-space homoge- 
neous solution. In the following, I will concentrate on the scalar perturbations alone (Si, R iy 
and will be neglected. 
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The vectors X M are now chosen to be the infinitesimal coordinate transformations re- 
quired to reduce the above general metric to some fixed form. For example, one choice (used 
by ABM) is 



X ll =[Q- 2HQ - B, 



d t Q] 



(6) 



This is exactly the infinitesimal coordinate transformation required to make the B and Q 
terms in the perturbation metric equal to zero, leaving only the diagonal 0, $ terms non- 
zero. But the equations of motion for these terms are also exactly the equations of motion 
for the gauge fixed metric, where the gauge fixed B and Q are set equal to zero. Although 
ABM mention a few such possible choices for the vectors X^, there are a large number which 
they do not mention because they restricted themselves to use only the local metric variables 
in defining X M , while the use of the matter variables and of non-local metric variables would 
allow an expanded set of possibilities. In fact, for any gauge fixing, one can find a vector 
X M which would implement that gauge fixing. For example, if we choose 



the gauge invariant scalar field perturbations, 6$, will be zero, as will the off diagonal tem- 
poral terms in the metric,!?. This choice of vector corresponds to the gauge fixing for the 
first order perturbations of the scalar field and scalar metric of Sip = 0, and 5g 0i = 0. The 
choice 



on the other hand puts the gauge invariant metric into the synchronous gauge. One problem 
with using such non-local gauge transformation is that because of the integration "constants" 
(actually functions of the spatial variables) in X, these choices still retain some residual gauge 
freedom, corresponding to those integration "constants" . 

The procedure which they follow is to transform their original general metric by means 
of these gauge transformations into what they call a gauge-invariant metric. However, this 
is simply the procedure of making an appropriate coordinate (gauge) transformation so as to 
put the metric into the appropriate gauge fixed form. They then calculate the second order 
fluctuations in this gauge fixed form, and obtain a result which they then claim to gauge 
invariant. It is, in the sense that the results have been calculated in a specific gauge, and 
thus that specific gauge allows no further gauge transformations, (or to put it in another 
language, the gauge transformation is undone by their reduction of the metric). 

The gauge invariant approach is thus identical to the gauge fixing approach and any 
problems or advantages of one are also problems and advantages of the other. 



To further the discussion, let me now solve the perturbation equations in the long wave- 
length limit. I will, with ABM, choose the longitudinal gauge fixing, which gives the following 
equations. The linear metric in this gauge is given by 




(7) 




(8) 



II. LONG WAVELENGTH SOLUTIONS 
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ds 2 = (1 + 2<&)dt 2 - a 2 (t)(l - 2<p)(dx 2 + dy 2 + dz 2 ) (9) 

where (fi(t,x,y, z),$(t,x,y, z) are the first order perturbation of the metric, and ip is the 
first order perturbation of the scalar field \& (called SPsi above). The first order Einstein's 
equation — Tij = for i ^ j gives 

didj{4> -<t>) = o, (10) 

which by rotational invariance gives 

(^,-iv%)(0-$) = O (11) 

for all i, j The solution is 

$ = + ^(t) x l + a(t)(:r 2 + y 2 + ^ 2 ) + (12) 

The term proportional to a corresponded to a perturbation toward one of the spherical or 
hyperbolic homogeneous space-times. The v and Kj terms can be removed by infinitesimal 
coordinate transformations. I will ignore these possibilities here. 
Einstein's equations are 

Efu, = G^ u — T^ u = (13) 
Selecting the two constraints, we have 

= 4 1} (0,^) = -6#0 -V4>- \v'i, - + ^V 2 (14) 



= El = ~di (2d t {a4>) + \ad t (y )i^ (15) 

The remaining (spatial diagonal) terms, are obtainable from these via the conservation 
equations E^ v . v — 0. 

In addition the linearized equation for the scalar field is 

P 1 = i> + 3Hijj - V 2 ip + V"(^ )*P - 4^o0 - QH^ (j) - 2^o0 = (16) 

In the above, \I/o(t) is the lowest order solution for the scalar field, a(t) is the lowest order 
scale factor, V(Psi) is the potential for the scalar field, and H = a/a is the background 
Hubble constant. These variables obey 

QH 2 = ^ 2 + V (17) 

U>1 = -2H (18) 

In the long wavelength limit, where we neglect the term V 2 0, the equation E\ t = and the 
linearized equation for the scalar field P 1 = are also linearized equations for homogeneous 
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perturbations of the universe. But at least two of those homogeneous solutions are simply 
coordinate transformations. These two correspond to multiplying the spatial variables by a 
constant (3 and the translating the time by an amount \{t). Ie, the new coordinates are 

t = t-X(t) (19) 

x i = (l + (3)x i (20) 

The metric now becomes 

ds 2 = (1 + 2X)dP - a(t) 2 (l -2(3 + 2HX)5 ij dx i dx j (21) 

Demanding that these transformations leave the metric in the longitudinal gauge with = $ 
we get 



which we can easily solve to give 



Thus, 



X = /3-HX (22) 



X(t) = Nadt + € (23) 



= t ( iladt+l\ f ^ = { _ H P_I^±1 + m ^ (24 ) 
\ a J a 

. ■ (3 J adt + e , 

$ = f(x l ) (25) 

(X 

must be solutions to E\ t = and P 1 = if we neglect the terms proportional to V 2 (ie in 
the long wavelength limit). The only remaining equation is the E\ = equation. But again, 
substituting our trial solution into that equation clearly satisfies it. 

Since the equation P 1 = is a second order temporal equation in ip, and the E tt = is 
a first temporal order equation in 0, we expect three linearly independent solutions to these 
two coupled linear equations. The above are two of them. However, solving E ti = for ip and 
substituting into E\ t = 0, we obtain a single second order equation in 0. Since the E\ t = 
equation uniquely determines -0 if is known, the linearized long wavelength equations have 
only two solutions, and we already know two of them. The above two solutions are the only 
solutions to the long wavelength equations. 

Note that without the E\ { = equation, there would have been three solutions. In the 
infinite wavelength limit, the E^ = equation is satisfied identically (because of the spatial 
derivative in that equation), and thus there exist three linearly independent solutions to 
the linearized equations. As I demonstrated above, two of them correspond to coordinate 
transformations of the original background solution, while the third would correspond to 
true a physical mode. It is thus interesting that both physical long wavelength solutions 
go to coordinate transformations in the infinite wavelength limit, while the actual physical 
infinite wavelength solution has no correspondence for large but finite wavelengths. 

Of course, having solved the problem in one gauge, one can determine the solution in any 
gauge. For interests sake, and because Grishchuk uses that gauge, let us look at the solutions 
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in synchronous gauge. Since we have the solutions in longitudinal gauge, the problem in 
synchronous gauge can be obtained by simply making the appropriate gauge transformation. 
The gauge transformation vector to go from longitudinal to synchronous gauge is given by 



(26) 



-J (f)(t,x,y,z)dt, di J 

The transformed metric is 

Sg(i)a = <*ff(i)ti = (27) 
Sg m = 2a 2 ((0 + H J <pdt)5^ + dtdj J l^j (28) 

Inserting the known solution for 0, i.e., 

# = /(,, V ,,)ft(i±££f!*) (29) 

we obtain 

5g (1)ij = 2a 2 (((3f(x, y, z) + k(x, y, z)H)5 ij + didjf(x, y, z) J e + ^ a dt dt ( 30 ) 

+didj(K(x,y,z) J a~ 2 dt + 'j(x,y,z))^ (31) 

k and 7 are temporal integration "constant" spatial functions. Both correspond to gauge 
transformations, illustrating the well known feature of the synchronous gauge that it does 
not completely specify the gauge. These terms could be removed by a gauge transformation 
which leaves the system in synchronous gauge. Note that the part of the metric which is 
the 4> term (whose contribution to the metric in synchronous gauge is proportional to 5ij) 
contains only the /3 parameter of the two parameter family of physical solutions. The other 
physical solution, given by the terms which depend on e, occurs only in the spatial derivative 
Q terms. Any second order equation for (/> would therefor pick up only one of the physical 
modes. Fortunately the solution which depends on e is one which dies out at long times 
in an expanding universe, and as a result their neglect would not be of importance in the 
late stages. Note that <fi is physically constant (up to a gauge transformation) for all times, 
agreeing with the contention of Grishchuk (where this term is given the name of h). 

Examining the long-wavelength perturbation in the scalar field, we find in synchronous 
gauge that 

i/> = * (t)K(x,y,z) (32) 

The only term remaining in the scalar field perturbation is a pure gauge term. Ie, syn- 
chronous gauge in the long wavelength limit is also a gauge in which the first order scalar 
field perturbations are zero, modulo a gauge transformation. 
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A. Grishchuk's Criticism 



Let me now comment on the controversy between Grishchuk and the rest of the com- 
munity. As has been emphasized by Grishchuk, who operates in synchronous gauge, the 
"growing mode" in synchronous gauge is constant for long wavelength modes. In the above, 
this is just the constant f3. He then argues that this indicates that there is no amplification 
of the scalar modes in inflation while the modes are outside the Hubble radius (ie long wave- 
length). Since in this regime the scalar field (matter) perturbations are zero, do not display 
any "amplification" either. His conclusion is that this indicates that the scalar modes are of 
the same size at the end of inflation (and in fact when they reenter the Hubble radius after 
reheating) as they were at earlier times when those modes left the Hubble radius (when the 
physical wavelength first became larger than the Hubble radius). Since the gravity wave 
modes also do not grow during inflation, this suggests to Grishchuk that both scalar and 
gravity wave modes should have roughly the same size at the end of inflation. This is in 
direct contrast with the standard lore, that in general the gravity wave modes are much 
smaller than the scalar modes. 

On the other hand, if one works in the longitudinal gauge, the growing mode of the 
gravitational metric parameter goes as (3H J adt/a — j3, which is (modulo decaying modes) 
identically zero if the expansion is exactly exponential, and is /3//i if the expansion goes as 
a power law a(t) oc t M . Furthermore the scalar field modes go as $(/3(t — to)/(/j, + 1)) f° r a 
power law increase in a. Since both \1/ and t increase as inflation continues, these scalar field 
modes grow from the start of inflation (where ^ and t — to are small) till the end (where \I/ 
is large.). 

But of course neither of these statements mean anything in themselves. It is clear that 
statements about the growth of perturbations are highly coordinate dependent statements, 
and thus contain no physics per se. The key point is that arguing whether or not a pertur- 
bation grows or does not grow is irrelevant if no statement is also made as to how large the 
perturbations were at the beginning of inflation. The whole argument between Grishchuk 
and others can be reduced to the contention on the one hand that it is ip in the longitudinal 
gauge which has a roughly vacuum quantum amplitude 1/vH (in units where h — 1) at 
the time when the modes cross the Hubble radius, while the other contends that it is the h 
mode, the diagonal term in the synchronous perturbation, 

ds 2 = dt 2 - (a 2 (l - 2h)5ij + Q 4:j )dx i dx j (33) 

which has that amplitude. (This view of the debate is weakened by noting that in ||, 
Grishchuk seems to adopt the same initial condition conventions as others). In short, it is 
crucial to decide how one will quantize the modes at the earliest times. By an appropriate 
choice any final result can be obtained. 

In also seems to an outsider that energy has been wasted about the details of the long 
wavelength calculations- Grishchuk contending the favourite "gauge independent variable" 
of the one group £ = + is zero in the long wavelength limit (it is not-it is just 

(5 in the above notation-see also [PHI , 0), while the others contend that he has not done 
his matching properly as the equation of state of the background changes (his h, as shown 
above can be chosen to be constant at all times, no matter what the potential V in scalar 
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field terms, or no matter what the equation of state in perfect fluid terms. It is however not 
clear that he has always chosen this gauge. ||) 

Trying to decide what the initial conditions are for the perturbations in the very early 
universe is a thorny issue. One tactic is to declare that these perturbations must be "vac- 
uum" perturbations. However the justification for this stance is somewhat weak, since by 
assumption (crucial for inflation to work at all) at least a part of the system is very far from 
its lowest energy state. The scalar field must have a very large value, and a large non-zero 
energy for inflation to proceed. However, to assume that it is far from equilibrium, while 
all other degrees of freedom at at their minimum energy is worrisome. But this is not the 
place to try to examine this issue in any detail. Thus I will use the assumption that the 
fluctuations are in some sort of minimum energy state. Energy is however a coordinate de- 
pendent quantity, and, because of the large non-equilibrium background field, there is a tight 
coupling between the gravitational degrees of freedom and the matter. One must therefor 
adopt a formalism which takes this into account. Fortunately, because of the background 
field, and the special coordinate system chosen for that field, one can define an energy for 
the fluctuations to lowest order in the those fluctuation. 



III. REDUCED HAMILTONIAN ACTION 

Let me therefor derive the reduced Hamiltonian for the scalar field perturbations in a 
flat FRW space-time (see also Mukhanov and Anderegg || and Garriga et. al Q). The 
Hamiltonian action for metric perturbations is 

S= [(7r ah 'y ab + vr^ - NH - N b H b )d x dt (34) 

where 

#0 = ^pir^ilaclcd - l/^ablcd) ~ yfiR + + 7^ A + V(V) (35) 

H a = Vl D b^= + K**,bJ ah (36) 



where 7 = det(j a b) and D a is the covariant derivative with respect to the spatial metric 7^. 
I will now write this action in terms of the first order perturbations for the scalar modes. 
Writing 

= + etp (37) 
7r* = a 3 ^ + $e7i> (38) 
7 ab = a 2 (l + 20)5 afe + g, a , 6 (39) 

n ab = -2Ha5 ab + e ((^^ + \^ 2 P Q W h + ^(-3V- 2 P Q - ^), c , d 5 ac 5 bd ^j (40) 

P^p, Pq have been chosen so as to be the true conjugate momenta to 0, Q. Ie, at second 
order 

/ ^ ab iab = e 2 J (P^ + P Q Q)d 3 xdt (41) 
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In addition, I choose 



N=l + ea ; N a = e/3 a (42) 

In the following, I will assume that all spatial dependence is of the form e ±tk ' x so V 2 becomes 
—k 2 . Substituting these into the action, and examining only terms at second order in e, one 
gets 

S = j + P Q q + P^ - aH 0{1) + PH^ >Q - H m d 3 xdt (43) 

where the () index indicates the order (in powers of e) of the term in question. Although 
I could write out the detailed form of the various terms, it would not be illuminating at 
this stage. To mimic the longitudinal gauge choice, I will choose my coordinates so that 
Q = Q = (3 = 0. From the variation equation of the action with respect to Pq, this results 
in the relation Pq = ^k 2 P^. Using the Hamiltonian equations for Pq andP^, derived by 
varying the action with respect to Q and (f>, one finds that this choice of gauge leads to the 
relation a = —<p. 

The resulting constraint equations are 

= if (i) = -i (~^o(a 3 H + 6a 3 HHP^)ij 
2H 

+8a 3 (3H 2 + 3H 2 H - Hk 2 )(f) + 2HHP 4> + 2^ #) (44) 
= H? lha = -k 2 Q^P^ + AHatjy - ai V>) (45) 

I then solved these for P^ and ip, and substituted the results back into the action. Clearly, 
the terms (the constraints) multiplying a and (3 vanish, and the action arises solely from 
the symplectic form and i?o(2) terms. The result 



r T 4k 2 1 / • 

s = / + ^ 2 ( HP * + {iSa3HH + 24a3H 



+ (288a 6 H 2 H + lUa 6 HH + lUk 2 a 4 H)(p 2 ) 



dt 



was still somewhat of a mess, so I defined the new variable w by 

= • (47) 

a 

Also, let I had to change define the momentum,P^, (which we note is no longer conjugate 
to 0) to a new one, P w , which was conjugate to w (ie, the momentum is conjugate if the 
action contains the term P w w with unit coefficient). The necessary transformation is 

p - = -^k ep * (48) 

The action finally reduces to 



S 



Ha „ H + H 2 ^ 2 2 / 2aH 2 H + a 2 HH + k 2 H\ 2 
P - w + IiP P ' + H W + 2HV { J 



dt (49) 
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Making another series of transformations 



™= (50) 



2kH 
P w = (Pw - w 



2H 2 



*£- (5!) 
-ifa 



and performing the appropriate temporal integrations by parts (in order to remove terms 
containing WPw or WW) we finally get 




-//a „ 2 / if 



if 



Note that this is just the Hamiltonian for a Harmonic oscillator with time dependent fre- 
quency. I could further reduce the action by changing time to conformal time r, and defining 
a new variable W = W/a. This would remove the \ dependence in the term proportional 
to k 2 , but the system is simple enough as it stands. 

This Hamiltonian action can be quantized in the usual manner. Assuming that at very 
early times the state for this variable is in the ground state for this Hamiltonian, one obtains 
that it is be the variable W which will have the vacuum quantum amplitude of fluctuations, 
which, using the WKB approximation at early times, gives \W\ ~ at those times. 

Furthermore one can solve the equation for w exactly if one makes certain assumptions 
about the background solution. Varying the action with respect to w and P w and eliminating 
P w gives an equation for w: 



if if \ k 2 
w-[H+^-2-\w + -w = (53) 



For long wavelengths, using our known solutions for 0, the solutions for w is — e — 
(3 J a(t)dt + j3a/H. (Substituting this into this equation for w verifies it as the general 
solution.) We now need to match this to the solutions for large k (or small time), using 
the now known amplitude for W at early times. The usual method is to match assuming 
that one is in De Sitter space-time at early times. However, I find this an uncomfortable 
procedure, as the equation for w, and the relation between w and W, the simple quantum 
variable, is singular as H goes to a constant, since both depend on H and H. Furthermore, 
both the equations for w and W are potentially singular in this limit. Instead I will assume 
that 

a(t) = t M , (54) 

with jj, a large constant. (The limit, /i — > oo, corresponds to De Sitter space.) Under this 
assumption for a, the equation for w becomes 

A* k 2 

w w H w = (55) 
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Defining the new variable r = — / * = , the conformal time, the equation for u> 

becomes 

d> + 2- — ^—d T w + k 2 w = (56) 
(/i - l)r 

which has as solutions Bessel functions 

„ = (*r)-*fc (-4 J ( 5 £±l y ,far) + B (57) 
The solution for W is then 

W = kJ^w = ky/ji((n - i)T)*tf=Vw (58) 
v a 

Now, as argued above, W will have amplitude of order ^a/2k at early times. This gives the 
equation for A, B of 

A = B = —=k^)(n-l)-^ (59) 

y/2fj, 



(A, B are actually quantum operators, and these expressions are shorthand for V < A 2 > 
and V< B 2 >.) Matching to the solution for long wavelengths (large t or small r), namely 

*> = -e- - J) = + ^TT)^ +1 ( 6 °) 

= A^% + B T -^ (61) 



Thus 



^ _ ,3 

K27T 



2^(^X1 + 1) 



fc2( M -l)(M=l)2T^T)( A(+ l)2( M 1 -l)r(I) 



(62) 



Since the term proportional to e dies out rapidly, I will not bother with giving its value. .This 
gives the amplitude for the quantum fluctuations during the period while the fluctuations are 
outside the Hubble radius. Note that they depend on ji. We can write ^fjl = J 1 . = 2^. 

The term in square brackets is a function of ji which is almost unity for large ji, but will give 
the slow change in the spectrum (k dependence) as \j/ changes slowly during the course of 
inflation. 

In conclusion, if one accepts the above procedure for determining the initial quantum 
fluctuations (an assumption which is not altogether beyond question) and if one accepts the 
validity of the linear theory for studying the determination of the initial amplitude for the 
fluctuations, then the answer which Grishchuk obtains for the final density fluctuations is 
wrong. 



11 



A. Comparison to other Reductions 



As mentioned, such reductions have been carried out previously. It may however be of 
interest to relate my original approach to this problem. Mukhanov, Feldman and Bran- 
denberger || report on such a reduction of the Lagrangian action. They derive a reduced 
Lagrangian action in terms of the longitudinal gauge variable 

v = aty - (*o/#» (63) 
Using the constraint, E^' = 0, this reduces to 

'UilA (64) 



a^ \H 

or, in terms of the reduced variable w, 

2H 



w (65) 
a\J-H 

My first reaction to this was to notice that in the long wavelength limit, w eliminates the long 
wavelength solution proportional to e and thus depends only on (3. Thus it would appear 
that the variable v does not capture the full set of physical solutions to the perturbation 
equations in the long wavelength limit, making it an unsuitable candidate for quantization. 
However, the relation is more subtle, and in particular it demonstrates another case where 
the k — > limit is not the same as the k = case. 



Examining equation 53, we notice that it can be written as a relation between w and v. 



In particular, we note that 



k 2 w = ^Ld t ^v (66) 



Thus, although the zeroth order e dependent solution for w does not contribute to v, the 
next k 2 order part of the solution does. There is thus a direct relation between the two 
linearly independent solutions for w and the two for v, despite the apparent contradiction 
when k = 0. 



IV. HIGHER ORDER PERTURBATIONS 

Let us now return to the question of whether or not the higher order long wavelength 
perturbations act to renormalise the cosmological constant, as discussed in ABM. Since to 
first order, the solutions locally (when spatial derivatives of the perturbations are neglected) 
look like simple coordinate transformations of the homogeneous solutions, their effect on 
the metric to higher order will again be that of simple coordinate transformations. Ie, in 
the long wavelength limit, the local effect of the perturbations will just be identical to a 
simple coordinate transformation of the background equations. As far as the local physics 
is concerned, the evolution of the universe is identical to that of a homogeneous universe, 
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with unrenormalised coupling constants. I.e., the effective cosmological constant will not be 
renormalized as far as local, sub-Hubble radius physics, is concerned. 

However, their analysis is not concerned with local physics. Rather what ABM argue 
is that, if we examine the average evolution of the universe in the large, these long wave- 
length perturbations act to alter the effective long wavelength evolution of the universe as 
a negative cosmological constant would. To examine this let us, again in the longitudinal 
gauge, examine the behaviour of the long wavelength perturbations. However, I will not 
follow their technique. In their approach, they wish to regard the first order perturbations 
as contributing an extra stress energy tensor to the zeroth order equations of motion (ie, 
the effective stress energy tensor of the first order equations renormalizes the zeroth order 
stress energy tensor.) I do not wish to follow this procedure as it raises delicate problems in 
consistency. Since by hypothesis, the zeroth order equations do not obey the zeroth order 
Einstein equations, a consistent derivation of the equations of motion becomes difficult. In- 
stead I will follow the procedure of consistently (though probably not convergent) expanding 
in a series of small perturbations. Ie, the metric will be assumed to have the form of 



where the equations of motion are then consistently expanded as a function of e. The 
effective stress energy tensor of the first order perturbations will act as the source for the 
second order perturbations in the metric and field. 

The biggest problem is that contribution to the second order terms from the first order 
perturbations (ie the effective energy momentum tensor for the second order perturbations) 
depend on the gauge chosen for the first order terms. Assume that we have two different 
gauge invariant formulations, determined by two separate vectors X M and X M . 

If in one system, the second order perturbations are given by Sg^)^u, then in the other 
gauge, the solutions will be (see appendix A) 



where A is the (first order) vector field which transforms the first gauge fixing to the second. 
A 2 is the second order part of this fixing, (note that in most of the work in ABM, was only 
concerned with the first order gauge fixing vectors. Since the second order equations were 
never solved, the need for the second order components of the X M was not there.) 

As I show in the appendix, by an appropriate choice of the first and second order gauge 
fixing, I can set the second order homogeneous perturbation in the metric to whatever value 
I wish. 

Let me be more definite. Let us first work in the ABM gauge, the longitudinal gauge. 

Going back to the observation that the first order solutions are, locally, just coordinate 
transformations, the reason for the apparent re-normalization of the cosmological constant 
is thus also clear. At each point the universe acts identically to the way it acts at each 
other point, except displaced in time by an amount which varies from place to place. If 
one now averages the universe over a fixed time slice, the averaged value of the expansion 
at fixed time will not be the same as the expansion rate at the averaged value of the time 
because of the non-linear nature of the expansion with time. Naively one would expect, 



9iw = 9oixv + egv„ + e g 2fll/ + ... 



(67) 
(68) 



&9(2)nv = Sg{2)tiv + 2£a&7(1)(/xi/ + ^A^aS'(O)^ + £a 2 5'(0)^ 



(69) 
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< a(r(t, x) >= a(< r(t, x) >) + a"(< r(t, x) >)(< r(t, x) 2 > — < r(t, x) > 2 ), where r is the 
uniform time (naively the time defined such that 5ip is zero along the r constant surfaces). 
The second term is thus present because of the non-linear relation between a and t. 

These conclusions about the also seem to be in agreement with the analysis of Salopek 
[IIJ who uses the Hamilton Jacobi methods to exactly solve for the evolution of the universe 
in the long wavelength limit, for restricted choices of the matter fields. His results seem to 
me also to suggest that the long wavelength evolution does not renormalise the cosmological 
constant. 

After completion of this work, the paper by Kodama and Hamazaki [Tr2[ was brought 
to my attention where they have found the solution to the long wavelength equations for 
a general multicomponent system. They also find that the long wavelength limit does not 
correspond to the homogeneous solutions for the scalar perturbations. 



APPENDIX A 

To derive the equation for the perturbations, let us first define what we mean by the 
expansion of the metric in the various orders. Consider the space of all solutions of Einstein's 
equations. This space can be defined by the set of tensor components g^ v {x) and fields 
4>(x). Now consider a path through this space of all solutions parameterized by an arbitrary 
parameter e. Thus we will have a set of functions g^vix, e), and we will choose the path and 
e such that at e = 0, the metric is the given background metric. For each value of e these 
are solutions of the equations, and thus we have 

G^u (jgap(x, e)) = T^(g al3 (x, e), <p(x, e)). (70) 

Since this equation is valid by assumption for all e, derivatives of these equations with 
respect to e will also be satisfied. Furthermore, assuming that these equations are analytic 
in e near e = 0, we can also create a power series expansion of the set of solutions along the 
path in e. 

Now consider another path through this space of solutions defined by an e dependent 
coordinate transformation 



x a = x a {x,e). (71) 



such that 



0x a dx^ 

g„u(x, e) = ^ dx J A^ e )> e ) ( 72 ) 

for some other path through the space of solutions defined by tildeg^ u (x, e). 
Now define 

_ dx a (x,e) dX^(x(x,e),e) 

where X^^x, e) is defined so that X fl (x(x, e), e) = x^ for all e. Also define the second order 
infinitesimal coordinate transformation by 



14 



^ = dri»(x,e) 
de 

Taking the second derivative of the equation [7^ with respect to e at e = 0, and assuming 
that x a (x, 0) = x a , we obtain the second order coordinate transformations. After some 
algebra, one finds that 

d 2 q uu (x,e = 0) dq uu (x,e = 0) „ „ . , „ dq uu (x,e = 0) „ „ . 

^ = ^ + e = 0) + 2£„ >- + C v Cr,g^{x, e = 0) 

(75) 

Ie, the second order metric components are not independent of the first order coordinate 
transformation. 

Let us consider the particular case where the first order solution = ^ e ' e ~°^ is 
evaluated in the longitudinal gauge (ie, is diagonal), and has spatial dependence with spatial 
wave- vector k. Assuming that k is very small (super-horizon modes), I will neglect terms 
proportional to k or k 2 . ( In the ABM paper, the first order k dependent terms would 
be zero because of the assumption that the fluctuations are statistically homogeneous, and 
have no preferred direction). Furthermore, let me assume that only n is non-zero, with 
arbitrary time dependence and with spatial dependence with wave vector k. Furthermore, 
let us look only at the effect of this first order gauge transformation on the second order 
spatially homogeneous terms. Then, 



and 



and 



C^g^v = 2rf5^5 0v + 2aar}°5 ij (76) 
C v (2h 00 + Cgoo) = r l °(2h 00 + 2ddotr] ) + A(h 00 + ?) )r) (77) 

C^kj + C v9ij ) = 2 V ° (fiij + rf(rf(aa + a 2 ) + ?) (aa))c%) (78) 

Using the gauge choice from the main paper, such that hij = 2a 2 $5jj, h 00 = 2phi, <fi = 
Phi ~ const, we have, neglecting terms of order k 2 , 

C v (2h 00 + C v g 00 ) = 2r/y + 4^° + 4(t) ) 2 (79) 
C^kj + C vgij ) = a 2 5 lj 2r ] °(2H(j) + (H + 2H 2 )r]°) + H(if) 2 (80) 

Note that all of the off diagonal terms due to this first order gauge transformation are zero 
because of the averaging over k. We know the gauge changes for the homogeneous modes, 
with C°(t) and Q = Xx l , which gives the changes 

£ C 3oo = 2(°C c g ij = (-2adC° + X)6 ij (81) 

The second order terms from the first order gauge transformations cannot be cast into this 
form. One thus has two independent functions of t, namely that produced by the first order 
inhomogeneous gauge transformation, and that caused by the second order homogeneous 
gauge transformation. By making appropriate choices of both the first order gauge with 
spatial wavenumber k and of the homogeneous second order gauge, and neglecting terms of 
order k 2 , one can set the two second order homogeneous terms in the metric (ie the 00 and 
the diagonal spatial components) to any value one desires, including 0. 



15 



ACKNOWLEDGEMENTS 



I would like to thank R. Brandenberger for discussion and for bringing the issue to my 
attention. I would also like to thank S. Mukhanov for extensive email in which he corrected 
my misconception about the relation between the MFB reduction and the one given in this 
paper. L. Grishchuk and I also have had extensive email discussions on the issues raised here. 
Finally I thank the Canadian Institute for Advanced Research and the Natural Sciences and 
Engineering Research Council for support during this work. The algebraic manipulations 
necessary in driving the equations and in reducing the Hamiltonian action were done with 
the help of GRTensorll under Maple V4. 



16 



REFERENCES 



[1] L.R.W. Abramo, R.H. Brandenberger, V.F. Mukhanov, Phys.Rev D56 (1997) 3248 
[2] S. Mukhanov, H. Feldman, R. Brandenberger, Phys Reports 215, 203 (1992) 
[3] S. Mukhanov, Anderegg Phys. Lett. B331, 30 (1994) 

[4] See also Jaume Garriga et al |astro-ph/ 9706229] "Canonical Quantization of Cosmo- 
logical Perturbations in the One Bubble Open Universe" for an application of similar 
Hamiltonian reduction to a hyperbolic spatial universe. 

[5] L. Grishchuk Phys Rev D50 7145 (1994) 

[6] L. Grishchuk gr-qc9801011| "Comment on the Influence of Cosmo logical Transitions on 

the Evolution of Density Perturbations." (1998) 
[7] J. Martin, D.J. Schwarz ^r-qc/9704049| "Influence of Cosmological Perturbations on the 

Evolution of Density Perturbations" (Phys Rev D to appear) 
[8] N. Deruelle,S. Mukhanov, Phys Rev D52 5549 (1995) 
[9] Robert Caldwell, Class.Quant.Grav.13 2437 (1996) 
[10] M. Goetz (to appear in Mon. Not. R. Aston. Soc) |astro-ph/9704271| " Validity of the 



'Conservation law' in the Evolution of Cosmological Perturbations" 
[11] For an alternative approach to the long wavelength perturbations using Hamilton Jacobi 
techniques see D.S. Salopek, J.R. Bond, Phys. Rev D42 3936 (1990); D.S. Salopek g 
| qc/ 9802025| "Hamilton- Jacobi Solutions for Strongly- Coupled Gravity and Matter" 
[12] Hideo Kodama, Takashi Hamazaki |gr-qc/9712045| " Evolution of Cosmological Perturba- 
tions in the Long Wavelength Limit." and Misao Sasaki, Takahiro Tanaka gr-qc/980117 
"Super Horizon Scale Dynamics of Multi-Scalar Inflation" 



17 



